Non-equilibrium sedimentation of colloids on the particle scale 
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Abstract 

We investigate sedimentation of model hard sphere-like colloidal dispersions confined in horizon- 
tal capillaries using laser scanning confocal microscopy, dynamical density functional theory, and 
Brownian dynamics computer simulations. For homogenized initial states we obtain quantitative 
agreement of the results from the respective approaches for the time evolution of the one-body 
density distribution and the osmotic pressure on the walls. We demonstrate that single particle in- 
formation can be obtained experimentally in systems that were initialized further out-of-equilibrium 
such that complex lateral patterns form. 
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The process of dispersed particles settling under gravity is a prevalent phenomenon in col- 
oidal dispersions, relevant for isolation of cells from blood [l| , separation of macromolecules 
B, industria. filtration S, and the shelf-.ife of paint Q . As demonstrated by the beautiful 
pioneering work of Perrin [5], particles with buoyant mass m that are exposed to a grav- 
itational field with acceleration g can be viewed as being subject to a (linearly) changing 
local chemical potential fi(z) = —mgz + const, where z is the height variable. This provides 
a means to access isotherms and e.g. study equilibrium phase behavior, provided that the 
gravitational length, £ g = k^T j{mg\ where ks is the Boltzmann constant, and T is tem- 
perature, is large compared to the inter-particle correlation length, but small compared to 
the vertical container size L . In equilibrium a pressure gradient that originates from 
an inhomogeneous density profile, p(z), balances gravity and for a dilute dispersion ideal 
gas-like behavior results, p(z) oc exp(— z/£ g ). In concentrated systems pronounced density 
oscillations, indicative of liquid layering, occur near the lower wall 8|. 

Understanding the time evolution of a given initial state towards sedimentation-diffusion 
equilibrium is a difficult problem due to the interplay of direct and (solvent-mediated) hy- 
drodynamic interactions between the particles jyj . An isolated sphere of radius a in a solvent 
of dynamical viscosity rj sediments at an average (Stokes) velocity Uq = Tomg, where the 
mobility is To = (67r?7a) _1 for no-slip boundary conditions on the sphere. Two timescales 
are of relevance, that of diffusion is the Brownian time it takes a particle to diffuse over 
a distance of one radius, tb = 3irrja 3 /(4ksT), and that of gravity is the time it takes a 
particle to sediment over one radius, ts = ct/(2m ), where a = 2a is the particle diameter. 
The ratio between both is the Peclet number, Pe = tb/ts = a/£ s , that characterizes the 
yp e of dynamics. Previous experiments were carried out using sedimenting granular matter 
where the effects of diffusion are small (Pe ^> 1), and using relatively small colloids 
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111 ] where diffusion dominates (Pe <C 1). Ste ady -state in the intermediate regime was 



recently investigated with computer simulations 12] showing that hydrodynamic backflow 
effects are relevant but depend little on Pe and that the velocity correlations are similar to 
those near equilibrium. 

In this Letter we investigate full non-equilibrium sedimentation of dispersions inside hori- 
zontal planar slits with plate distance L and Pe ~ 1 using experiment, theory and simulation, 
all of which are capable of resolving the microscopic structure of the dispersion on the par- 
ticle level. We obtain quantitative agreement of results from the respective approaches for 
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the strongly inhomogeneous density profiles. This contributes to our understanding of the 
effects of particle correlations and confinement on the non-equilibrium time evolution of sed- 
imentation and could be relevant for transport of cells in microfluidic devices bacteria 
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in groundwater [17| and in porous media 

Our experimental system 19] consists of micron-sized, sterically stabilised, fluorescently 
labelled Poly methyl methacrylate particles [^(J (PMMA), dispersed in a mixture of cyclo- 
hexyl bromide (CHB) and cis-decalin with 260 pM tetrabutyl ammonium bromide (TBAB) 
salt to screen the (weak) electrostatic interactions between the colloids, resulting in a Debye 
screening length of ~100 nm [3]; see Tab. [I] for a summary of all relevant parameters. We 
tailored the density difference between solvent and particles to control the buoyant mass m 
and hence Pe. Initially dry capillaries were filled with the suspension; this is a very rapid 
process due to large capillary forces between the walls and the solvent. The capillaries were 
laid flat such that the colloids sediment towards the bottom wall, and imaged with a confocal 
laser scanning microscope, see Fig. UK for an illustration. Confocal scans were made in three 
dimensions and time series are presented as slices in either vertical (xz, see Fig. [Tb-d') or 
horizontal (xy) planes. Coordinates of the particles were tracked in three dimensions [21] , 
yielding the colloid density profile, p(z,t), as a function of the altitude z and time t. We 
present results for L/a = 12 (system A) and L/a = 38.5 (system B). 

Experimental results for the time evolution of the sedimentation profile for system A 
are shown in Fig. [2^l. The inital density distribution is almost flat, save for a small peak 
at the lower wall. This peak grows in time and two more layers appear. The final profile 
exhibits clear second and third layers and decays smoothly for higher altitudes. We have 
observed no pronounced changes in the lateral (xy) structure, nor strong fluctuations in 
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and contrasting with the case of 



particle velocity, in accord with the findings o : 
strongly driven (steady-state) sedimentation jlO ]. 

A recent and very promising theory for the time evolution of the one-body density, p(r, t), 
where r is the position coordinate, is dynamical density functional theory [22j, [23( (DDFT): 

where T is the mobility, F[p] is taken to be the equilibrium Helmholtz free energy functional, 
F[p] = ksT f drp(r) [ln(p(r, t)A 3 ) — 1] + F exc [p] + f drp(r)V ext (r), where the first term is the 
Helmholtz free energy of the ideal gas, A is the (irrelevant) thermal wavelength, -F exc [p] is 
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the excess (over ideal gas) contribution to the free energy, and the external potential in 
the present case is V ex t( r ) = ksTz/^g for a/2 < z < L — a/2 and oo otherwise. The 
interparticle interactions are described by the excess free energy functional F exc [p], which 
we take to be that of hard spheres of diameter a c g > a to model the (strongly screened) 
electrostatic repulsion. We use Rosenfeld's successful fundamental-measures theory [24| to 
approximate F exc [p\. As a simple treatment of hydrodynamic interactions we take into 
account the reduction of T upon increasing density with the result by Hayakawa and Ichiki 



25), r ffl (0)/r o = (1 - 0) 3 /[l + 20 + 1.4920(1 - 0) 3 ], where = 4vrpa 3 /3 is the colloid 



packing fraction, and we generalize to strongly inhomogeneous situations by using the local 
packing fraction 0(r,t), obtained as the convolution of the bare density profile, p(r,t), 



with a weight function characteristic of the particle volume [24j. We determine values 
for C, g and a e s such that the theoretical result for the equilibrium sedimentation profile 
[obtained from solving SF[p]/5p(r) = p for p(r), where p is the chemical potential] matches 
its experimental counterpart. The resulting agreement is very good indeed, see the top pair 
of curves in Fig. Eti- In addition we have checked that values for a c s are consistent with 
our experimental results for the radial distribution function (not shown) and that the values 
for C, g are compatible with those obtained from estimating m from sample composition. 
We assume p(r,t) = p(z,t), take p(z,t = 0) to be the equilibrium profile for hard spheres 
between hard walls without gravity, and integrate Eq. ([1]) forward in time. The results for 
p(z, t) compare very well with the experimental data, including the time evolution of the 
layering at the bottom which was not captured in earlier treatments 9], where monotonic 
density profiles were obtained - as is appropriate on a macroscopic scale. 

We further test the DDFT by comparing results for p(z, t) to those from BD simulations. 
To correct for the finite resolution of the simulations, we have used an intrinsic clock to 
measure time: The clock ticks only when a trial move is accepted, but not when a trial 
move is rejected due to hard core overlap. We have extensively tested this approach against 
time step extrapolation [26[] and are confident that it is accurate using time steps of 1CT 3 tb, 
hence offering significant advantages in computational efficiency We use a constant mobility 
Thi(0o) both in BD and DDFT, with O being the overall packing fraction. The agreement 
between results from theory and simulation is remarkable, see Fig. [2b, establishing for the 
first time that DDFT is accurate for describing Brownian dynamics of hard spheres in 
three dimensions; this is both of fundamental and practical importance due to the high 
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computational efficiency of the scheme. 

System B serves to investigate the crossover from strongly confined to macroscopic sys- 
tems. In the early stages we find a large region at intermediate heights where the density 
remains nearly constant, reminiscent of batch settling [2^], as well as a pronounced upper 
interface of the sediment towards an essentially particle- free supernatant, see Fig. [2b. For 
long times pronounced layering at the bottom of the container develops; we have analyzed 
the lateral structure in the experiment and have ascertained that these layers remain liquid- 
like, with no signs of crystallization, which is expected to occur for yet more layers 
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achieve very good agreement of the experimental and theoretical results for p(z,t), albeit 
with a small mismatch in time scales: The theoretical results are shown for a value of tb 
that is 20% larger than the experimental value tb = 29.3 s, obtained from measurements 
of particle size and viscosity. The discrepancy might be due to particle size polydispersity 
(~ 5%) becoming more important in the case of the larger value of Pe in system B, and not 
taken into account in the DDFT. 

The detailed structural information contained in p(z, t) enables us to investigate the time 
evolution of important macroscopic quantities such as the osmotic pressure [8| exerted by 
the particles on the bottom and top walls, Ph{t) = ksTp(z = cr/2,t) and Pt{t) = ksTp(z = 
L — a/2,t), respectively. Results shown in Fig. [2bl indicate a sharp initial increase in Pb{t) 
with a crossover to slow relaxation for long times towards its equilibrium value, whereas 
Pt{t) very rapidly decays (not shown). Note that in equilibrium P^ — P t = ksTu/^g, where 
v is the number density per area. 

In order to drive the system further out-of-equilibrium we have altered the initial state by 
turning equilibrated samples upside-down, such that the colloids preferentially reside at the 
top of the capillary. For this system £ 9 fh 0.35cr and L fx 18a, such that Pe ~ 1.6. Despite the 
fact that these parameters are comparable to those considered above, sedimentation proceeds 
in an entirely different fashion. A strong degree of lateral inhomogeneity developsjte9J] , see 
Fig. El In the horizontal (xy) plane, see Fig. [3b, a network-like structure emerges |30j. The 
structure formation is reminiscent of that resulting from the Rayleigh- Taylor instabililty in 
two layered liquids, given that the steep initial density gradient resembles that of a liquid- 
liquid interface. There is likely a relationship between this phenomenon and the swirls 
observed in sedimenting granular matter [10]. A theoretical investigation of this instability 
requires an explicit treatment of the lateral x and y coordinates and a more sophisticated 
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handling of the solvent hydrodynamics [31] than only through a density-dependent mobility; 
these tasks are beyond the scope of the present work. We believe that DDFT, when extended 
along these lines, could be a very promising candidate for studying such phenomena, both 
on the grounds of the high accuracy in simpler situations (demonstrated above) and that it 
also has been used successfully in the past to describe pattern formation in processes such 
as spinodal decomposition [32 ]. 

In conclusion we have presented the first measurements on the single-particle level of 
sedimentation in model colloidal dispersions. Unlike previous macroscopic descriptions that 
apply when the gravitational length is much larger than the particle size |9|, we have consid- 
ered cases where both lengths are similar in magnitude. We have demonstrated that DDFT 
gives accurate results as compared to computer simulation results and that it describes 
laterally homogenous sedimentation very well as compared to experimental data. Besides 
systems that exhibit lateral pattern formation it would be interesting to investigate in future 



work larger systems where crystallization occurs as well c 
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O va 2 a L £ g /a L/a t b (s) a cS /a Pe 


A 
B 


0.15 3.624 3.5 50 0.8 12(+1) 57.9 1.05 0.625 
0.142 10.46 2.8 100 0.45 38.5(+l) 29.3 1.1 1.11 



TABLE I: Overview of parameters for systems A and B: Total hard sphere packing fraction <f) , 
scaled density per lateral area va 2 , where 4>o = ira 3 v/(6L), particle diameter a (in microns) 
as obtained from static light scattering, capillary spacing L (in microns), scaled gravitational 
length £, g /a, scaled plate separation distance L/a, Brownian time tb (in seconds), scaled effective 
(electrostatic) hard sphere diameter a e g/a, and Peclet number Pe. The values for £ g /a and a e s/a 
are obtained by fitting theoretical results for equilibrium density profiles of hard spheres to those 
from experiment. The lateral size of the capillaries was 1 mm x 50 mm. 
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FIG. 1: Non-equilibrium sedimentation of colloids, (a) Illustration of colloidal spheres of diameter 
a under the influence of gravity g and vertically confined between two walls with separation distance 
L. The vertical coordinate is z; the horizontal coordinates are x and y. (b), (c), (d), Time series 
of confocal micrographs of system A (see Tab. |T] for parameters) taken in the (vertical) xz plane 
at times t/rs = 3,26,200, respectively. The scale bars denote 20 ^m; the horizontal lines indicate 
the position of the walls. 
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FIG. 2: Time evolution of the density distribution, (a) The local packing fraction cf)(z, t) = 
p(z, t)ira 3 /6 as a function of the (scaled) height coordinate z/a for (scaled) times t/rs = 
3, 7, 14, 39, 100 (from bottom to top) for system A. Shown are results from experiment (full lines) 
and from DDFT using a density-dependent mobility (dashed lines) . The curves are vertically offset 
by 0.2 for clarity, (b) Same as (a) but obtained from BD computer simulations (full lines) and 
DDFT (dashed lines); in both approaches a constant mobility is used, (c) Same as (a) but for 
system B for times t/rs = 7,26,44,63,255 (from bottom to top). Time t = is set such that 
the results for the earliest time match; the theoretical curves are scaled with a value of tb that is 
20% larger than the experimental value. The experimental result for t = 255t# is compared to the 
theoretical equilibrium profile, (d) The osmotic pressure that the particles exert on the lower wall, 
Pb(t), as a function of the scaled time, t/rs, as obtained from DDFT for system A (lower solid 
line) and system B (upper solid line). The asymptotic values for t — > oo are indicated by dotted 
lines. Results are obtained using a density-dependent mobility. 
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FIG. 3: Rayleigh- Taylor-like instability. An equilibrated sample is turned upside down and left to 
sediment ((a), vertical xz plane at time t = 2tb). A strong fingering-like inhomogeneity develops 
quickly ((b), vertical xz plane at t = 5tb), accompanied by maze-like lateral pattern formation 
((c) horizonal xy plane at time t = 9tb). Subsequently particle diffusion equilibrates the system 
((d) vertical xz-plane at time t = lira), similar to the time evolution shown in Fig. [Tb-d. The 
dashedjine in (b) indicates the height of the xy plane in (c). The scale bars denote 20 ^m; see 



also 



3C|. 
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